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Abstract. 

An expression for a third-order link integral of three magnetic fields is 
presented. It is a topological invariant and therefore an invariant of ideal 
magnetohydrodynamics. The integral generalizes existing expressions for third- 
order invariants which are obtained from the Massey triple product, where the 
three fields are restricted to isolated flux tubes. The derivation and interpretation 
of the invariant shows a close relationship with the well-known magnetic hclicity, 
which is a second-order topological invariant. Using gauge fields with an SU(2) 
symmetry, helicity and the new third-order invariant originate from the same 
identity, an identity which relates the second Chern class and the Chern-Simons 
three-form. We present an explicit example of three magnetic fields with non- 
disjunct support. These fields, derived from a vacuum Yang-Mills field with 
a non-vanishing winding number, possess a third-order linkage detected by our 
invariant. 



1. Introduction 

The topological structure of magnetic fields is an important subject in plasma physics. 
There, among other issues, it is related to the problem of stability of a plasma and to its 
energy content. Fields with an enormous wealth of entangled, braided or knotted held 
lines exist for example in the solar atmosphere. Note that the topological complexity 
of these solar magnetic fields is only revealed, if first, one takes into account that 
the observed loops anchored in the photosphere are closed by subphotospheric fields, 
and second, that already simple toroidal equilibria contain many different knotted 
and linked field lines. The simplest examples are the so called torus knots, which are 
formed by field lines where the quotient of the number of windings around the core 
and the torus axis is rational. 

In general, magnetic fields in plasmas are not static, but evolve due to the 
motion of the plasma. The evolution of solar as well as most astrophysical 
magnetic fields is given in good approximation by the induction equation of ideal 
magnetohydrodynamics (IMHD) 

^B-Vx (vxB) = 0, (I) 
at 

which shows that the field can be considered as frozen-in with respect to the plasma 
velocity v. The approximation of IMHD which leads to this law is valid as long as the 
evolution does not lead to small scale structures, e.g. thin current sheets. 
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The ideal induction equation guarantees the conservation of the topology of field 
lines under the flow of v, i.e. every linkage or knottedness of magnetic flux is preserved. 
Mathematically speaking, the flow of v is a differentiable isotopy of the space under 
consideration. It maps the field lines of B at time to to a topologically equivalent set of 
field lines for any later time t > to. Let us note that, following the usual terminology 
in plasma physics, the term 'topological equivalent' is used here in the sense of a 
diffeomorphic isotopy. 

In order to describe the structure of magnetic fields, it is desirable to have 
measures of complexity. These measures should be topological, i.e. they should be 
invariant under an arbitrary isotopy of the magnetic field, and therefore invariant 
under the ideal induction equation. An example of a topological measure for magnetic 
fields is the magnetic helicity, a quantity which has attracted a great deal of attention 
in recent years (see e.g. Brown et al 1999). Magnetic helicity, which measures the 
linkage of magnetic flux, is only a lowest order topological measure. It fails for instance 
to detect the interlocking of magnetic flux tubes in form of the Borromean rings or 
the Whitehead link (see Figure |l|) . The total magnetic helicity of these configurations 
vanishes just as it does for three or two unlinked flux tubes. Both configurations must 
possess a higher order linkage or knottedness of magnetic flux which is not detected 
by magnetic helicity. This naturally raises the question whether corresponding higher 
order measures similar to magnetic helicity exist which are sensitive to these linkages. 
Here we would like to remark that the configurations shown in Figure |l| are highly 
idealized. In any real plasma we would not find this pure linkage but a mixture of 
different types of linkages, each of which is to be measured by a different integral. 

In knot theory different invariants are known which distinguish e.g. the 
Borromean rings or the Whitehead link from unlinked rings. The problem is that 
invariants used in physical applications, e.g. in magnetohydrodynamics, should be 
expressed in terms of observable quantities, in our setting in terms of the magnetic 
field B. Up to now only the helicity, which is related to the Gaufi linking number, 
has been formulated as an invariant for magnetic fields in a satisfactory manner. 
As was recognized first by Monastyrsky and Sasorov (1987) and independently by 
Berger (1990) and Evans and Berger (1992), the link invariants based on so-called 
higher Massey products (see Massey 1958, 1969, Kraines 1966, Fenn 1983) can be 
written as invariants applicable to magnetic flux tubes. Similar to helicity, they only 
involve the magnetic fields and can be expressed as volume integrals over the space in 
consideration. Their disadvantage is that their usage is restricted to magnetic fields 
confined to isolated flux tubes. In addition, these flux tubes must not possess a linkage 




Figure 1. The Borromean rings (left) and the Whitehead link (right). 
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lower than the linkage which is measured. 

In this paper we present a generalized third-order invariant for three magnetic 
fields not confined to isolated flux tubes. In the case of isolated flux tubes this invariant 
coincides with the invariant known from the Massey triple product. Using gauge 
fields in the context of an SU (2) symmetry, the generalized invariant can be shown to 
originate from the same equation as helicity. Therefore, we will first recapitulate some 
basic facts about magnetic helicity before we turn to our main subject, the third-order 
link invariant. 



2. Magnetic helicity 

Magnetic helicity of a field B with arbitrary vector potential A is defined as 

H{B)= [ A-Bd 3 i, B-n\ dv = 0, (2) 

Jv 

which can readily be shown to be gauge invariant if no magnetic flux crosses the 
boundary of the volume. Since it is quadratic in magnetic flux it is often referred to 
as a second-order topological invariant. Magnetic helicity measures the total mutual 
linkage of magnetic flux. This interpretation can be motivated if we envisage a 
simple system of two isolated and closed flux tubes U\ and U 2 with vanishingly small 
cross-section. The latter condition ensures that the (self-)helicities of the flux tubes 
vanish. Moffatt (1969) has shown that for this configuration the helicity integral yields 
H(B) =24> 1 4> 2 lk(U!,U 2 ), where lk(U u U 2 ) is the Gaufi linking number (Gaufi 1867) 
of the two tubes and where <f>i is the magnetic flux in the tube Ui. Introducing an 
asymptotic linking number, Arnol'd (1974) was able to extend this interpretation to 
the generic case where field lines are not necessarily closed (see also Arnol'd & Khesin 
1998, Part III, §4). 

Similar to helicity we can introduce the more general cross-helicity of two 
magnetic fields Bi and B2. For a simply connected volume V and provided that 
Bi • n\gy = B2 • n\gy =0we define 

J7(Bi,B 2 ):= / A 1 -B 2 d 3 s= / Bi-A 2 d 3 :r. (3) 
Jv Jv 

The boundary conditions ensure that this is again a gauge invariant quantity. To 
see that the two integrals on the right-hand side are equivalent, note that they differ 
only by a surface integral J(Ai x A 2 ) ■ nd 2 x. This can be shown to vanish using 
the equivalence of B; • n\gy = with A, x n\gv — in a certain gauge and for a 



simply connected volume V, as proved in Appendix A. Since both integrals are gauge 



invariant this proves the equality for any gauge. 

H (Bi, B 2 ) measures purely the cross linkage of flux among the two fields. Applied 
to our system of two isolated closed flux tubes with fields Bi + B 2 = B in the 
corresponding tubes this invariant yields 

H(B±,3 2 ) = 2 cf>i<f>2 lk(£/i, U 2 ), (4) 

which is now valid without the assumption of vanishingly small tube cross-sections. 

The significance of magnetic helicity arises from the fact that it is invariant in 
IMHD. Using merely the homogeneous Maxwell equations we obtain 

d t (A • B) + V ■ (A Q B + E x A) = -2E • B, (5) 
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which describes the time evolution of helicity density A • B. Here A is the electric 
potential and E the electric field. The term A TS + E x A is to be interpreted as 
a helicity current, and — 2E • B as a source term. In an ideal plasma, i.e. with 
E + v x B = 0, the source term vanishes and the helicity current can be written as 

AoB + ExA = A B + v(A-B)-B(A-v). (6) 

Therefore equation (^|) takes the form: 

9 t (A-B) + V-(v(A-B))=V-(xB), (7) 

with x = — + A • v. Elsasser (1956) already noticed that a particular gauge can be 
found such that x = 0. Using either this gauge or an arbitrary gauge together with 
the boundary condition B • n| a y = this last equation implies the conservation of 
helicity in a comoving volume for an ideal plasma, since 

J A • B d 3 x = J d t (A • B) + V • (v (A ■ B)) d 3 x = J %B • nd 2 x = 0. (8) 

v v av 

The invariance of integral (|j) was first stated by Woltjer (1958). 



3. A third-order invariant from the Chern-Simons three-form 

In this section we construct a third-order invariant which, under conditions specified 
below, yields an invariant for three magnetic fields. The derivation is based on some 
basic knowledge in differential geometry found e.g. in Frankel (1997). 

We have noted before that equation (||) can be derived purely from the 
homogeneous Maxwell equations. Written in differential forms it reads 

d{A A dA) = F A F, (9) 

where A is the one-form potential of the field F. The right-hand side of this equation 
is one of two (pseudo-) scalar Lorentz invariants that can be constructed from the 
field tensor. We can interpret this equation as a special case of a general result in the 
theory of Chern forms, namely the exactness of the second Chern form 

2 

dtr(A AdA— iq-A A A A A) = tr(F A F). (10) 
3 

In this equation A and F are a matrix valued one- and two-form. To be more precise 
they take their values in the Lie-Algebra g of the structure group. In Yang-Mills 
theory this is the symmetry group of the interaction under consideration, with coupling 
constant q. On the vector space q the trace defines a natural scalar product. For to, 
9 £ Q, tr(w A 6) = Uij A 6ji, where the indices denote matrix components. 

Equation © holds for an arbitrary, not necessarily Abelian, field strength 
F = dA — iqA A A. The three-form on the left-hand side is known as the Chern- 
Simons three-form. For the case of electrodynamics, i.e. for the Abelian structure 
group U(l), F is given by F — dA, since AAA vanishes and equation ( [To|) reduces 
to (Q) . In the non- Abelian case equation ([Io|) splits into a real and imaginary part 

tr(dAAdA) =dtr(AAdA), (11) 

2 tr(A A A A dA) = dtr(-A A ^4 A A), (12) 
3 

As we will see in the following a third-order invariant can be derived from identity 
( p^|) for the special case of the structure group SU(2). 
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Working with an SU(2) structure group it is appropriate to choose the Pauli 
matrices a i: i=l,2,3, as a basis for the Lie Algebra su(2). All results, however, are 
independent of this choice. The gauge potential A and field strength F now have three 
components 

A = a l A\ F = a i F i , (13) 

where the summation convention over repeated indices is assumed. Let us note that 
in the following we will refer to these fields as Yang-Mills fields, although they do not 
necessarily satisfy the Yang-Mills equation. Using the identities for Pauli matrices 

<jj(Tk = i ejki 01 + Sjkl and aj o^ai = iejkil, (14) 

equations (|ll]) and ( |l2| ) read 

dA i A dA l = d(A i A dA*), (15) 

e ijk dA* A A> A A k = i d{e l]k A l A A* A A k ). (16) 

If we now interpret the three components of the Yang-Mills potential A as three 
independent potentials of three electromagnetic fields F l EM = dA l , the first identity 
states the helicity conservation (in IMHD) for the sum of the self-helicities of the three 
individual fields F EM , similar to the electrodynamic case. The second identity is new. 
For convenience we introduce the two-form G 1 = A 2 A A 3 , here on Ikf 4 , and cyclic 
permutations of it. Then we can write the second identity as 

~d(A i AG i ) = F% M AG i . (17) 

To complete the analogy of this equation with equation (^|) we have to rewrite it in 
the language of three- vectors. Therefore we represent the one-form A 1 by the time 
component A l and the three-vector A 1 of the corresponding four-vector. The two-form 
G 1 is identified with the vector pair (G E , G B ) = (A 3 , A 2 — AqA 3 , A 2 x A 3 ), equivalent 
to the identification of Fem with the three- vector pair (E,-B). Cyclic permutations 
immediately lead to corresponding pairs for G 2 and G 3 . Using these conventions, the 
left- and right-hand side of equation (|l7j ) read respectively 

id(A< A CP) = -1 (d t (A* • G^) + V • (A^ + A* x G^)) vol 4 
= -(a t (A 1 -G 1 3 ) + V.(^G 4 B ))vol 4 

and 

Fem A G l = (E l -G l B -B l - G E )vo\ 4 = (E 4 • G^ - A l V • G^)vol 4 . 

Thus, identity ( p7| ) is equivalent to 

^(A 1 • G^) + V • (4JG' B ) = -E 4 • G l B + A l V ■ G B , (18) 

which shows a similar structure as equation (||) in the case of helicity. It describes the 
time evolution of the density A 1 • A 2 x A 3 with its current A l G B and source term 
— W ■ G' B + A l V ■ G l B . This is the basis for the following theorem. 

Theorem: Let B 1 , B 2 and B 3 be three magnetic fields with potentials A 1 satisfying 
E'+vxB' = 0. The integral over a volume V C 1R 3 

ff( 3 )(B 1 ,B 2 ,B 3 ) := [ A 1 -A 2 xA 3 d 3 a; (19) 
Jv 

is a gauge invariant, conserved quantity, if 
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(i) the potentials obey V • (A 4 x A?) — for all i,j = 1, 2, 3, 

(n) i/ie potentials obey the boundary condition A 4 x n|gy = / or n being the normal 
vector to the boundary of the integration volume V . 

Proof: Let us first remark that condition @ of the theorem implies B 4 • n\gv = 
as shown in Appendix A. It is therefore consistent with con dition (fl) sin ce = 
(B 4 • A 3 — A 4 ■ W)\q V = V • (A* x A J )|ay. Moreover, we show in Appendix A that for 
a simply connected volume with B 4 ■ n\gv = condition @ can always be satisfied. 
In order to prove the invariance of 7J (3) we observe that in an ideal dynamics equation 
(Q~q) can be written as 



^(A 1 • A 2 x A 3 



V 
V 



(v(A* • A 2 x A 3 )) 



iv(A* 



Gb) 



(A 4 — A 4 ■ v)V ■ G 



B - 



where the last term vanishes due to condition (|) . 

Integrating over the volume V yields the total time derivative of on the left- 
hand side, while the right-hand side can be converted into a surface integral, analogous 
to equation (||) in the case of helicity. The surface integral vanishes since condition 
([]) of the theorem implies = on the boundary of t he volume V. This, together 
with the gauge invariance of shown in Appendix B, completes the proof. 



4. Interpretation of the invariant 

It is interesting to note that the new third-order invariant comes on an equal footing as 
the conservation of helicity, since both invariants where derived from the same identity 
([To|). However, contrary to the conservation of helicity, the third-order integral cannot 
be applied to a single magnetic field, but requires a triplet of fields. Thus we have to 
interpret this integral in the sense of the cross-helicity rather than the total helicity. 
A forthcoming paper will deal with the question of how a single magnetic field might 
be split into a triplet of fields with the required properties V • G l = 0, thereby linking 
the given cross third-order invariant and a total third-order invariant. 

There is another way of looking at the new third-order invariant. By writing 

H^(B\B 2 ,B 3 ) = [ A 1 • G l B d 3 x = H(B 1 , G#) , (20) 
Jv 

the integral is to be interpreted as the cross-helicity of the two divergence- free fields B 1 
and G]j. Note that the boundary conditions for the cross-helicity, namely B 1 • n = 
and Gg ■ n = 0, are fulfilled due to condition (||) of the theorem. From this new 
interpretation the condition V • Gg = is an obvious requirement analogous to 
V • B — 0. Furthermore, the symmetry of leads to 

i/(B 1 ,Gi,)=i/(B 2 ,G 2 3 )=i/(B 3 ,G 3 3 ), (21) 

which reveals the additional conditions V • G 2 B = and V ■ G B = 0. Let us note that 
this interpretation does not simplify the direct calculation of the third-order invariant. 
It is still necessary to determine the fields G^ which are not independent of the chosen 
representatives A 4 . 

Third-order linking integrals for magnetic fields have been constructed from the 
Massey triple product already by Monastyrsky and Sasorov (1987), Berger (1990) and 
also Ruzmaikin and Akhmetiev (1994). However, these constructions are limited to 
cases where the three fields are confined to three isolated and mutually unlinked flux 
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tubes with disjunct support. It is in fact easy to see that for this special case their 
invariants coincide with the integral ( |l9| ) given above. An explicit proof is given in 



Appendix C . In particular, it is worth noting that for fields with mutually disjunct 



support the condition V • G l B = implies that the cross-helicity of all pairs of fields 
vanishes, i.e. their flux tubes have to be mutually unlinked. 

For a set of three arbitrary magnetic fields V ■ = cannot always be satisfied. 
To show that there are examples for which this can be satisfied, beyond the cases of 
three fields with mutually disjunct support, we give an explicit example in the next 
section. 



5. Example of three magnetic fields with a third-order linkage 

In this section we want to give an example of three magnetic fields not confined to 
flux tubes, which firstly allow for one- form potentials that obey d(A l A i 3 ) = and 
where secondly the integral invariant ( |l9|) yields a non-trivial result. The existence of 
such an example proves that the new invariant is indeed a generalization of the third- 
order invariant derived from a Massey triple product which was applicable merely to 
unlinked flux tubes. The fields we construct show an extraordinary high symmetry. 
For this reason they are interesting in their own right. 

The idea to construct three fields A 1 on M 4 which obey d(A l A A 7 ) = comes 
from Yang-Mills theory: An SU(2) Yang-Mills field 

F = dA— iqA A A (22) 

can, in view of equation (113) and the identities for Pauli matrices (|l4|), be written as 



F i = dA 1 + q e ijk A j A A k . (23) 

By taking the exterior derivative of F, 

dF* =qe ijk d(A j f\ A k ), (24) 

we immediately observe that dF — is a sufficient condition for all d(A l A A?) to 
vanish. In the special case of a vacuum Yang-Mills field, i.e. F = 0, the requirement 
dF = is trivially fulfilled. If we now reinterpret the three components of the Yang- 
Mills potential as potentials of three independent magnetic fields, we have constructed 
an example field configuration to which the invariant (|l9| ) can be applied. 

5.1. Yang- Mills potentials of a vacuum field with non-vanishing winding number 

An SU (2) Yang-Mills vacuum field is now constructed on a time slice M 3 of M A using 
the mapping g : Ft 3 — > SU(2) (see e.g. Frankel 1997, Itzykson and Zuber 1980) 



S(x) = exp 



ITT x J a 



(25) 



Interpreted as a gauge transformation of an SU(2) classical vacuum, i.e. with 
vanishing connection lu = 0, g(x) gives rise to the pure gauge connection 

^ = .g- 1 (x)d g (x), (26) 

and the Yang-Mills potential one-form A = ^uj reads 

A*<Ti = V'Md^x). (27) 
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At this point we want to remark that the vacuum winding number of to, which is 
defined to be the degree of the map <?, is W(g) — 1. An important consequence of this 
non-trivial winding number will be a non-trivial value of the invariant H^ 3 \ as can be 
seen in equations ([32]) and ( |33| ) below. 

In order to explicitly calculate the one-form potential A given by the last 
expression we use that g(x), as an element of SU(2), has the form 

ff(x) = exp (^nA = cos (-) - i^n? sin (-) , (28) 



where n = /n and where n is a unit vector in 1R 3 with coordinate components 
fiP = ic J '/||x||. A comparison with equation ( |25| ) shows 

27r||x|| 

/ = ; ■ 29 



Substituting equation (g8|) into (g7j) we obtain after some calculation 

A i = — [n l d/ + dn 4 sin / + e y fe n J dn fc (cos / - 1)] . (30) 
2q 

5. 2. Magnetic fields constructed from Yang-Mills potentials 

As mentioned above we now interpret the three components of the non-Abelian 
pure gauge Yang-Mills potential A as potentials of three independent, non- vanishing 
magnetic fields. It is sufficient to consider only one of the three potentials A\ since 
due to n % = a; l /||x|| and the cyclic symmetry of equation ( |30| ) in the indices k 
all three fields can be obtained from just one field by rotations that map the x'-axes 
on one another. From equation ([50|) we calculate the vector potential A 3 in spherical 
coordinates x — r sin i? cos ip, y = r sin ■& sin ip, and z — r cos i?. Using unit basis vectors 
e r , e&, e. v and fixing a value for the "coupling constant" of q = 2, we find 

. , — 27TCOS?? siwd ( 2irr \ 2 sin?? 2 / irr \ 

■ e r H sm sin e v . (31) 



" (l + r2)3/2 r r \VTT^J r ^ 

The corresponding magnetic field can be calculated from B 3 = V x A 3 . It is easy to 
check that the fields B z are well defined and scale as ||B J (x)|| — » An 2 for r — > and 
as ||B l (x)|| — > r -4 for r — > oo. Hence, they have no singularity and decay sufficiently 
fast for large radii. 

Let us note that the fields B 1 are highly symmetric and similar. Looking at the 
vector potential A 3 we observe that it is independent of the variable <p, therefore B 3 
is invariant under rotations leaving the x 3 -axis fixed. Since the potentials ( |30| ) are 
cyclic in i,j, k it follows that each B* field is invariant under rotations about the x l 
axis. We have pointed out before that a rotation that maps the Euclidian basis vector 
field ei to e2 also maps B 1 to B 2 etc. Furthermore, B 3 is similar to the total field 
B = B 1 +B 2 +B 3 in the following sense: Let R be a rotation that maps the Euclidian 
basis vector e 3 to the vector VSi?^) = ei + e 2 + e 3 , then B = V3 i?(B 3 ). 

The magnetic fields B J are only of interest for us if their third-order invariant does 
not vanish. Explicitly calculating A 1 for i — 1, 2, 3 we find, using the main Theorem, 

ff <3><B',B*,B3) . jf" HTpfrr) - * - < 32 » 

The fact that this integral is non-vanishing proves, that the constructed invariant 
cannot only be applied to all cases for which we where able to calculate the already 
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existing invariant, i.e. to three mutually unlinked flux tubes, but also to examples of 
triples of fields not having disjunct support, ft is thus a true generalization of the 
existing invariant known from the Massey triple product. 

As we have pointed out before, H 1 --^ is related to the vacuum winding number 
W(g) of the connection ui — g~ 1 &g. We easily find (see also Frankel 1997) 

W[g) := ^ J tr{g- l dg A g^dg A g^dg) 

l — [ e l]k A l A A j A A k (33) 



96tt 2 



v 



= -^ (3)(Bl ' B2 ' B3) ' 
where the trace term is usually referred to as the Cartan three-form on SU(2). 

In the general case the cross-helicities of three magnetic fields, for which we are 
able to find potentials such that V • (A 1 xA J ) = B l A- 7 — A* • B J = 0, do not have to 
vanish. In our example we can easily verify that they do vanish, i.e. H(B l , B J ) = for 
i j. Of more interest are the three non-trivial sclf-hclicities. If a triple of magnetic 
fields is derived from a Yang-Mills vacuum, equation ( |24| ) together with F = implies 

dA* + ~ ei jk Ai A A k = 0. 
2 

Using the definitions B % = dA 1 and G % — A> A A k we find for (i, j, k) cyclic 

B l = -A 1 A A k = -G\ 
Thus, for i — 1,2,3, we observe that 

H(B\B l )= [ A*-B l d 3 x = - [ A ,; • G'd 3 i = -ff (3) (B\ B 2 , B 3 ). 
Jv Jv 

Therefore the self-helicities are equal to the value of the third-order invariant. This is 

a peculiarity of all magnetic field triples derived from an SU (2) Yang-Mills vacuum. 

In our analysis of the three example magnetic fields we now turn our attention to 

the topological structure of the fields and the linkage of individual field lines. Figure || 

and U show numerically integrated field lines, where the starting points for integration 

are indicated by the foot points of the arrows that give the field line direction. We 

observe that all field lines are closed and have an elliptical shape. Figures | and | 

visualize the toroidal structure of the individual fields at the example of B 3 . Using 

A 3 = A 3 r e r + A^ + A 3 e v and ^A 3 = it follows that B 3 can be written 

B 3 = V x A 3 = -J-V(r sinM 3 ) x e v + B%e v . (34) 
rsmv v v 

Therefore, the field lines of B 3 lie on (^-invariant toroidal surfaces described by 

rsiniM 3 = const. Figure || shows the poloidal B 3 -field and contour lines for three 

different values of rsinzM 3 . A toroidal surface with rsin$A 3 = — 1.8 and four field 

lines on it is drawn in Figure |3[ The central field line, sitting within all tori is 

characterized by £? 3 = B# = 0. In view of the last equation, this is equivalent to 

V(r shifM 3 ) = which yields r = ^/l/3 and ■& = n/2. We observe that all field 

lines wind around the central field line exactly once. From this and the toroidal 

structure of B 3 we can conclude that any two arbitrary field lines I and I' of B 3 have 

a Gaufi linkage lk(l, V) = 1. 

Finally let us discuss the linking properties among field lines of different fields. 

To give an example, one field line of each field is plotted in Figure ^[ The symmetric 
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appearance is due to the choice of symmetric starting points for the field line 
integration. As was stressed above, the magnetic fields B l can be obtained from 
one another by cyclic permutations of the Cartesian coordinates x % . In the same way 
the integration starting points for the field lines 7 1 , j 2 and 7 3 where chosen to be the 
cyclic permuted coordinate triples (0, 0.8, 0), (0, 0, 0.8) and (0.8, 0, 0). It is interesting 
that the total linkage of the set of field lines is lk( 7 1 , 7 2 ) + lk( 7 2 , 7 3 ) + lk(7 3 , 7 1 ) = -3. 
Even though we have seen that the mutual cross-helicities of all three fields vanish, 




x' 



Figure 2. The vector plot shows the 
projection of the B 3 -field onto the x 1 - 
z 3 -plane. The arrow lengths are not to 
scale. The solid lines are contour levels 
with values —0.5, —1, —1.8 of the plotted 
density distribution r sin 1?/1 3 . 



Figure 3. The drawn torus is the 
surface at which rsin$A 3 = —1.8. All 
field lines of B 3 lie on toroidal surfaces 
described by r sin $A 3 = const. Four 
field lines lying on this torus are drawn. 
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their individual field lines, in general, are linked pairwise. To be more precise: If 
we e.g. fix one field line 7 3 of B 3 , then all field lines of B 1 and B 2 are either linked 
with 7 3 exactly once or they intersect 7 3 twice. For reasons of symmetry, there exists 
for each field line 7 1 of B 1 a field line 7 1 of B 1 , such that we find the total linkage 
lk(7 1 ,7 3 ) +lk(7 1 ,7 3 ) = 0. Hence, the cross-helicity ^(B 1 ^ 3 ) vanishes, which in 
turn implies i^B^B 3 ) = 0. Figure || shows such field lines 7 1 , 7 1 and 7 3 with 
(0,0.8,0), (0,-0.8,0) and (0.8,0,0) as their respective starting points for the field 
line integration. Finally let us remark that in the same way as we obtained 7 1 we 
can obtain field lines 7 2 and 7 s , here with integration starting points (0,-0.8,0), 
(0,0,-0.8) and (—0.8,0,0). Together, these three field lines yield a configuration 
complementary to the one shown in Figure |], which now has a total linkage of +3. 

6. Conclusions 

An integral expression has been presented which generalizes the third-order invariant 
known from the Massey triple product, to an invariant not limited to mutually unlinked 
flux tubes, if the involved fields allow for potentials that obey V ■ (A 1 x A 3 ) = 
for i 7^ j. An example shows that the new invariant is a true generalization. 
In our derivation helicity and emerge from the same general identity, which 

involves the Chern-Simons three-form in the context of an SU (2) gauge symmetry. 
Whether this identity leads to further results for other gauge groups has not yet been 
investigated, but it is clear that only expressions quadratic and cubic in magnetic flux 
can be obtained. The constructed invariant is to be seen as a "cross-linkage" of three 
fields. It still remains to clarify whether or how a total third-order invariant can be 
constructed and whether this is possible with the help of a cross third-order linkage 
such as in the case of helicity. There might e.g. exist a subdivision of a single field 
into three components such that the total third-order linkage is determined by the 
cross-linkage alone. Unfortunately, the antisymmetry of iJ^^B 1 , B 2 , B 3 ) seems to be 
one of the key problems for a further generalization analogous to helicity. 
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Appendix A. Equivalence of boundary conditions 

We prove for a simply connected volume V the equivalence of the boundary conditions 
B • n\gv = and A x n\gv = 0. 

First note that A x n\gy = implies B • n|gy = 0: Locally on the boundary we 
can write n = V/3, where /3(x) is a scalar function defined such that /3(x) = (3q defines 
the boundary dV. Thus A x n\gy = implies A\gv = aV/3 for some a(x). Then 
= Va x V/3 and therefore B • n\g V = 0. 

To prove the reverse we start with an arbitrary vector potential A which 
will in general have a non- vanishing component A|| tangential to the surface dV. 
We can express A|| as a one-form a defined only on dV. Then the assumption 
B ■ n\gv = (V x A||) • n = written in differential forms reads da = on dV. 
From V being simply connected it follows that dV has the same homotopy type as 
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the two-sphere S 2 . But since the cohomology vector space iJ 1 (S' 2 ; M) — 0, all closed 
one-forms are exact. Therefore there exists a scalar function ip on dV such that 
a = dtp. This in turn implies that a gauge exists such that An|gy = and thus 
A x n| a y = 0. 

Appendix B. Gauge invariance of the third-order invariant 

We now prove that the integral invariant ([l9]) is unchanged under all gauge 
transformations A 1 — > A 1 ' — A 1 + V0 l for i = 1,2,3, which obey the following 
two conditions: First we require 

V • G 1 = V • G 1 ' = 0, (B.l) 

where as before we define G' = A- 7 x A k and G 1 ' = A? x A fc ' for cyclic indices k. 
Second, the gauge transformations must respect the boundary condition 

A* x n\ dv = A 1 ' x n\ dv = 0, (B.2) 

where n is a normal vector to the boundary dV of the integration volume V. 

It is easily checked that a general gauge transformation that leaves the condition 
V • G* = unchanged for i = 1,2,3 has to be a simultaneous gauge transformation of 
all three fields. Substituting A 1 + V</> 1 for A 1 our invariant of equation ([l9]) changes 
according to 

H (3) _> H (3)' — f (A 1 _j_ V 0l) . [( A 2 + V ^2) x ( A 3 + v ^3jj d 3 x ( B .3) 

Jv 

= f A 1 - (A 2 x A 3 )d 3 a; 
Jv 

+ [ [V0 1 • (A 2 x A 3 ) + \7<f> 2 • (A 3 x A 1 ) + V0 3 • (A 1 x A 2 )] d 3 x 
Jv 

+ [ [A 1 - (V0 2 x X7(f> 3 ) + A 2 - (V</> 3 x V0 1 ) + A 3 - (V^ 1 x V</» 2 )] d 3 a; 
Jv 

+ [ V^ 1 - {\7<j> 2 x \7cj) 3 )d 3 x. 
Jv 

We have to show that for gauge transformations respecting equations (BT) and (B.2), 
HW =H&, i.e. the sum of the last three integrals in equation flB.4| ) has to vanish. 
We can rewrite the first integral as 

[ V</> 4 • G* d 3 a; = [ <j) l G l • ndi - / <j) l (V ■ G l ) A 3 x = 0, 

Jv JdV Jv 

which vanishes since A 1 x n\gy = implies G* • n\gv = 0. To show that the second 
integral vanishes we use the identities 

V(j) 2 x V0 3 = V x O 2 V0 3 ) = -V x {(j) 3 V(j) 2 ) (B.4) 
= -V x {<t> 2 V4? - <t?V<t> 2 ) (B.5) 

and 

V • G 1 ' - V ■ G 1 = V0 3 ■ B 2 - V0 2 • B 3 , 
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as well as expressions obtained by cyclic permutations of the indices (1,2,3). 
Substituting the first identity into the second integral we obtain 



A 1 • (-V x (0 2 V0 3 - 3 V^ 2 )) + (1, 2, 3) cyclic 



d 3 x 



(A 1 x (< 

(v x a 1 ; 



2 v</> 3 



J V(?!> 2 )) • n+ (1,2,3) cyclic 
- 3 V</> 2 ) + (1,2,3) cyclic 



d x 



d 3 x 



6 1 (B 3 V</. 2 - B 2 V0 3 ) + (1,2,3) cyclic 
6 X (V ■ G 1 ' - V • G 1 ) + (1, 2, 3) cyclic 



d A x 



d 3 x 



= 0. 



The surface i nteg ral vanishes due to condition (BJ2) and the last volume integral due 
to conditi on ( B.l). Finally we can see that the third integral vanishes by substituting 
equation (BJj) for the term V</> 2 x Vcj) 3 

V0 1 • (V</> 2 x V^ 3 ) d 3 x = i / (V0 1 ■ V x (<?(> 2 V0 3 - 3 V<?!> 2 )) d 3 x 



2 V0 3 



3 V</> 2 ) x V0 1 • nd 2 x 



In the last step we used that due to condition (B^) the gradients V0 1 , V0 2 and V</) 3 
have to be parallel to n. This completes the proof. 



Appendix C. Equivalence of the link integrals for disjunct flux tubes 

The equivalence of the third-order link integrals as given by Monastyrsky and Sasorov 
(1987), Berger (1990) and Ruzmaikin and Akhmetiev (1994) for three disjunct 
and mutually unlinked flux tubes Ui with the integral (^) is shown as follows. 
Monastyrsky and Sasorov gave an integral which corresponds to the Massey triple 
product and reads in vector notation: 

/ (AixFi-A 3 xF 3 )d 2 izi, with V x Fi = Gi := A 2 x A 3 , for 1,2,3 cyclic. (C.l) 

The integration is taken over the surface of tube U\ . Cyclic permutations of indices 



in (C.l) yield equivalent expressions. Note that the G^ in this representation are 
evaluated only outside the tubes Ui where V • G^ = for any gauge. To convert this 
integral into a volume integral over the whole space, one has to evaluate Gi and F^ 
on Ui and therefore encounters the problem of V • (Aj x A^)^ ^ for an arbitrary 
gauge. To overcome this problem Berger defined G^ within the flux tubes as 

!Aj x A k - '!> , / H, on U j: 
Aj x A fe + ^yBfc on Uk, with i, j, k cyclic and V$(j)j = A j \ Ui .(C.2) 
Aj x A fc else, 

The resulting volume integral 

B 1 • (F a - $ (1)2 A 3 ) d 3 x (C.3) 

v 
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is equivalent to (C.l), as shown in Berger (1990). Using the same construction with 
potentials Ruzmaikin and Akhmetiev (1994) have rewritten ( p.3[ ) in a more 

symmetric form. 

Now, instead of using the additional potentials we can just as well use the 

special gauge 

A, -> Ai := Aj - V$ (j)i - V$ (fc)l , 

which implies Ai\i/ j — f or i ^ j. In other words, the corresponding new are 
set to zero and definition (|C.2| ) implie s con dition (i) of the Theorem in section \ 
Furthermore, suppressing tildes, ( C.3 ) turns into 

/ B 1 xF 1 d 3 a;= / A x x G x d 3 x + / (A 1 xF 1 )-nd 2 x 

JV JV JdV 



Ai x Gi d A x + (n x Ai) • Fi d x 

\ JdV 

A l x Gi d 3 x , (C.4) 



which finally shows the equivalence of the integrals flC.l|) and (|19| ) for the case of three 
pairwise unlinked flux tubes. 



References 

Arnol'd V I 1974 Proc. Summer School in Differential Equations (Erevan) Armenian SSR Acad. Sci. 

[English translation: 1986 Sel. Math. Sov. 5 327-45] 
Arnol'd V I and Khesin B A 1998 Topological Methods in Hydrodynamics Applied Mathematical 

Sciences vol 125 (New York: Springer- Verlag) 
Berger M A 1990 Third-order link integrals J. Phys. A: Math. Gen. 23 2787-93 

Brown M R, Canfield R C and Pevtsov A A (eds) 1999 Magnetic Helicity in Space and Laboratory 

Plasmas Geophysical Monographs vol 111 (Washington: American Geophysical Union) 
Elsasser W M 1956 Reviews of Modern Physics 28 135 

Evans N W and Berger M A 1992 A hierarchy of linking integrals Topological aspects of fluids and 
plasmas Nato ASI Series E vol 218 ed H K Moffatt et al (Dordrecht: Kluwer Academic Publisher) 
pp 237-48 

Fenn R A 1983 Techniques of geometric topology London Mathematical Society Lecture Note Series 

vol 57 (Cambridge: Cambridge University Press) 
Frankel T 1997 The Geometry of Physics. An Introduction. (Cambridge: Cambridge University Press) 
Gaufi C F 1867 Werke vol 5 (Gottingen: Konigliche Gesellschaft der Wissenschaften) p 602 
Itzykson C and Zuber J-B 1980 Quantum Field Theory. (New York: McGraw-Hill) 
Kraines D 1966 Massey higher products Trans. Am. Math. Soc. 124 431-49 

Massey W S 1958 Some higher order cohomology operations Symp. Int. Topologia Algebraica, Mexico 
(UNESCO) pp 145-54 

1969 Higher order linking numbers Conf. on Algebraic Topology, Univ. Illinois at Chicago Circle, 

June 1968 ed V Gugenheim pp 174-205 

Reprinted in: 1998 J. of Knot Theory and Its Ram. 7 No.3 393-414. 
Moffatt H K 1969 Journal of Fluid Mechanics 35 117-29 

Monastyrsky M I and Sasorov P V 1987 Topological invariants in magnetohydrodynamics Sov. Phys. 
JETP 66 (4) 683-688 

Ruzmaikin A and Akhmetiev P 1994 Topological invariants of magnetic fields, and the effect of 

reconnection Phys. Plasmas 1 331-336 
Woltjer L 1958 Proc. Nat. Acad. Sci. 44 489 



